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Abstract. WRA 751 is an evolved massive star in our Galaxy 
closely resembling Luminous Blue Variable stars (LBVs). It is 
surrounded by a nitrogen enriched nebula of about 23 " diame- 
ter. A comparative study of the nebula's morphology and kine- 
matics is presented, it supports — together with spectroscopi- 
cal evidence — the classification of WRA 751 as a LBV. Images 
show that the nebula consists of a nearly spherical shell as well 
as a bipolar-like structure north and south of its main body, the 
Northern and Southern Caps. 

In contrast to the almost spherical appearance of the main 
body of the nebula, the kinematics shows a deviation even of 
this part from a classical spherical expansion pattern. From the 
present data it can be concluded that the main body expands 
asymmetrically (central expansion velocity ~ 26 km s _1 ), with 
a thicker shell at the back side. A bump-like structure can be 
found to the west of the central star. In addition to the main 
body, bipolar kinematic components can be identified with the 
morphologically classified Caps. These results put WRA 751 
into the class of LBVs which are surrounded by a nebula with 
bipolar components, albeit considerably less pronounced than, 
for instance, in the classical bipolar LBVs r] Car and HR Car. 

Key words: Stars: evolution - Stars: individual: WRA 751 - 
Stars: mass-loss - ISM: bubbles: jets and outflows 



1. Introduction 

With today's initial mass function, stars with masses of M 



50 — 100 Mq and luminosities of L 



10 



5-6 



Lq top the 



Hertzsprung-Russell diagram (HRD) and represent the most 
massive stars known. These massive stars show a remark- 
able evolutionary behavior (e.g., Schaller et al. 1992, Stothers 
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& Chin 1996). After spending a 'normal' life as O stars 
on the main sequence they evolve towards cooler tempera- 
tures while entering a phase with very high mass loss (up 
to 1O _4 M0 yr _1 ), they become Luminous Blue Variables 
(LBVs). This phase starts when the stars reach the Humphreys- 
Davidson limit (Humphreys & Davidson 1979, 1994, Langer 
1994) in the HRD. Analyzing HRDs of the Galaxy and the 
LMC, Humphreys (1978, 1979) and Humphreys & Davidson 
(1979) found a lack of very luminous red supergiants. Apper- 
ently the most massive stars do not evolve into red supergiants 
but instead their evolution is reversed towards the blue super- 
giant part in the HRD. The turning points form the empirical 
Humphreys-Davidson limit. Around this Humphreys-Davidson 
limit in the HRD LBVs are found. One of the most promi- 
nent characteristics of the unstable LBV phase are strong stellar 
winds and possible giant eruptions, which lead to the peeling 
off of parts of the stellar envelope and the formation of small 
circumstellar nebulae, the so-called LBV nebulae (LBVN, e.g. 
Notaet al. 1995). 

Only very few LBVs are known in our Galaxy (8 classified 
and candidate objects) and a few in other galaxies, for instance 
in the LMC, SMC, M31, and M33 (Humphreys & Davidson 
1994). Altogether nearly 40 LBVs and good candidates are cur- 
rently known. Not all of them, but many show circumstellar 
nebulae. In this paper we analyze the LBVN around the galac- 
tic star WRA 751 and put it into context with the well known 
and better studied LBVs r] Carinae and HR Carinae. 

WRA 751 (= He 3-591) was first identified as a possible 
Wolf-Rayet star by Henize (see Roberts 1962) and then ap- 
peared in the Carlson and Henize (1979) sample of southern 
peculiar emission-line stars, where the authors already recog- 
nized the strong [Fe II] lines and classified it as Bep. In addi- 
tion they proposed a similarity between WRA 751 and HR Car 
(which in the meantime is known to be a LBV). Using spectro- 
scopic and photometric observations Hu et al. (1990) noted for 
the first time that WRA 751 is a good LBV candidate. In addi- 
tion to the characteristically strong [Fe II] lines this star shows 
an irregular light variation typical for LBVs. In their analysis 
the star was classified as 09.5 with Mboi = — 9.6 m . They de- 
rived T eff = 30 000 K, E(B-V) = 1.8, M = 50 M , and a 
distance of about 5 kpc. 



2 



K. Weis: The asymmetric nebula around WRA751 



One year later Hutsemekers & Van Drom (1991a) found 
that WRA 751 was surrounded by a nebula roughly 22" in di- 
ameter. They found the nebula to be nearly spherically sym- 
metric and to show a diffuse structure. Their low-dispersion 
spectra revealed a nebula of low excitation with an electron 
temperature of T e?ne b < 15 000 K, an electron density of n e ~ 
400 cm -3 , and an expansion velocity of v exp = 26 km s _1 . As- 
suming a distance of 7 kpc (derived from radial velocities and 
making use of the galactic rotation curve) their radius of the 
nebula is 0.38 pc. 

Analysis of infrared data (see de Winter 1992) show a NIR 
excess from which an estimate for the mass loss of M ^ 
10 -5.5...-6 Mg) yr -i with a velocity of v wind = 500 kms -1 
were derived. They found a cool dusty circumstellar shell with 
strong emission in the FIR to surround WRA 751. In 1992, 
from interstellar and circumstellar reddening, van Genderen 
measured a distance of 4 — 5 kpc for WRA 75 1 as a lower limit. 
Garcia-Lario et al. (1998) re-determined the characteristic pa- 
rameters and found T e ^ ~ 25 000 K, v exp = 24kms _1 , and 
n e - 200 cm- 3 . 

In summary, the published observations and their interpre- 
tations point towards WRA 75 1 being a LBV star with a slowly 
expanding nebula. 

2. Observation and data reduction 

2.1. Imaging 

To compare kinematics and morphology we used images of the 
nebula around WRA 751 taken with the Super Seeing Imager 
(SUSI) on the New Technology Telescope (NTT) of the Euro- 
pean Southern Observatory (ESO). The images were retrieved 
from the ESO NTT archive and reduced in a standard way, us- 
ing flat-fields and bias frames from the same observing night. 
Besides of the emission-line images in H a and [N II] a broad 
band red continuum filter image was used for the continuum 
subtraction in order to be able to distinguish the pure line emis- 
sion from the continuum)]. Fig. |l| shows a 1050 s exposure H a 
image (not continuum subtracted) with a round field of view, 
about 50" in diameter. From the image, we determined the see- 
ing as 0"8. As the image was taken using a coronograph, an 
occulting bar obscures the central star and some parts of the 
nebula. A non-occulted image of the nebula can be found in 
Hutsemekers & van Drom (1991a). 

2.2. Long- slit echelle spectroscopy 

To analyze the nebula around WRA 751 in more detail we 
preformed high-resolution long-slit echelle spectroscopy. We 
used the echelle spectrograph on the 4 m telescope at the Cerro 
Tololo Inter-American Observatory. Inserting a post-slit H a fil- 
ter (6563/75 A) and replacing the cross-disperser by a flat mir- 
ror we selected the H a line and two [N II] lines at 6548 A and 

1 Due to missing documentation in the NTT archive and especially 
the lack of information in the header of the images a more detailed 
specification of the filters can not be given here. 
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Fig. 1. ESO NTT archive image of the nebula around 
WRA 751, taken with SUSI using an H a filter and a corono- 
graphic mask (see also Nota 1998). The diameter of the field of 
view is about 50". The Northern and Southern Caps are indi- 
cated. A north-east vector marks the celestial orientation. 

6583 A. We picked the 791 mm echelle grating and a slit- 
width of 150 /i (corresponding to 1"), which lead to an instru- 
mental FWHM at the U a line of about 10 km s _1 . 

All data were recorded with the long focus red camera and 
a 2048 x 2048 Tek2 CCD.The pixel size was 0.08 Apixel" 1 
along the dispersion and / /26 pixel -1 on the spatial axis. Vi- 
gnetting limited the slit length to ~ 4' . Seeing was ~ 2" 
during the observations and the weather was not photometric. 
Thorium- Argon comparison lamp frames were taken for wave- 
length calibration and geometric distortion correction. 

The slit was oriented in two different position angles (PA) 
perpendicular to each other (Fig. ^). We observed 5 positions 
with PA = 140°. One position was centered on the star (in 
our nomenclature Slit center), two positions were offset to the 
north by 8" and 12" from the central star (Slit 8N and Slit 12N) 
and two positions were offset by 8" and 16" south of the cen- 
tral star (Slit 8S and Slit 16S). For PA = 50° we observed at 
two positions offset from the star by 9" to the north and south, 
respectively (Slit 9N and Slit 9S). The left column of Fig. | 
shows the echellograms of the slits at PA = 140°, Fig. [3] those 
at PA = 50°. All echellograms cover a spectral range of 75 A 
centered on H a and extend about 2' in the spatial direction, 
centered on the projected position of the central star. The top 
of each echellogram in Fig. [3] points to the north-west, and in 
Fig. [j]to the south-west. Some telluric lines are visible and were 
used to improve the absolute wavelength calibration. 
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3. Morphology of the nebula around WRA 751 

The first image taken of the nebula around WRA 751 was pub- 
lished by Hutsemekers & van Drom (1991a). They discovered 
a spherical diffuse nebula about 22" in size. We retrieved and 
re-analyzed [N II] and H a images retrieved from the ESO NTT 
archive in order to compare the morphology and the kinematics 
of the nebula. Figure jj] shows that the nebula around WRA 751 
indeed appears nearly round and almost spherically symmet- 
ric (see also Nota 1998). We measure a diameter of the main 
body of the nebula of 22'/ 8, which corresponds to 0.50 pc as- 
suming the lower limit distance by van Genderen et al. (1992; 
~ 4.5 kpc). Appearing rather homogeneous in its surface den- 
sity and closely attached to the central star, it resembles very 
much a Stromgren sphere. However, this is in contrast to the 
morphological appearance of most other LBV nebulae like that 
of AG Car or HR Car (e.g. Hutsemekers & van Drom 1991b, 
Nota et al. 1995, Smith et al. 1997) which show a more intense 
detached ring structure, or as in the case of HR Car at least 
parts of a bipolar shell (Weis et al. 1997, Nota et al. 1997). 
This is mainly due to an increased matter density along the line 
of sight when looking through the edge of a shell. 

A closer inspection of the images allows us to locate spa- 
tial variations in the surface density of the nebula around 
WRA 75 1 : a brighter circular half- shell shows up in the eastern 
part, while — in comparison — in the western part the brightness 
decreases by about a factor of two. Even more remarkable are 
two extensions of the spherical shell, one nearly exactly to the 
north at a position angle of 340° which we will call the North- 
ern Cap (Fig. [l]) and the other to the south (Southern Cap) at 
PA = 165°. This means that the Caps are almost along an axis 
through the star. Morphologically both appear roughly trian- 
gular in shape. The surface brightness of the Northern Cap is 
somewhat larger, that of the Southern Cap somewhat smaller 
than the low surface part (the west side) of the central spheri- 
cal nebula. The Northern Cap extends by 2'! 6 beyond the main 
body of the nebula, the Southern Cap by 4 / /7. 

4. Kinematic indication of a bipolar structure 

Using high-resolution long-slit echelle spectra an analysis of 
the kinematics of the LBVN around WRA 751 was performed. 
With the 7 slit positions our mapping covers the entire neb- 
ula. The right columns of both, Figs. || and ^ show the 
position- velocity diagrams (pv-diagrams) corresponding to the 
echellograms in the same line. For the pv-diagrams we used 
the stronger [Nil] line at 65 83 A which is less contaminated 
by background emission. Additionally all measurements were 
compared to those of H a to ensure that no kinematic differ- 
ence exists between the two lines. All velocities were corrected 
to the Local Standard of Rest (LSR) system. The 0" position 
in all pv -diagrams corresponds to the projected position of the 
central star (WRA 751) onto the slit, from there positive offsets 
are to the north-west for the PA = 140° slits (or south-west 
for PA = 50°) and negative to the south-east (or north-east, 
respectively). All spectra (except Slit 16S) clearly resolve the 




Fig. 2. Same image as in Fig. [I] with the slit positions overlay ed: 
two slits (9N and 9S) are oriented at PA = 50°, all others at 
PA = 140°. 

Doppler ellipse, and indicate a predominantly spherical expan- 
sion of the nebula. The largest radial expansion velocity v exp 
was found in the central position (0") of Slit 9S (PA = 50°) 
with a value of 29.4kms _1 ± 2kms _1 , the second largest in 
Slit 8N at 28.3kms _1 ± 2kms _1 . A similar expansion ve- 
locity was derived at the same position in the central slit at 
26.3kms -1 ± 2kms _1 . Since this slit crosses the star in the 
center, the velocity measurements are less certain at the 0" po- 
sition. There v exp might be even higher as can be estimated 
from interpolating for the missing data points in the pv-diagram 
of Fig. ||. If a nebula is spherically expanding, the largest ex- 
pansion velocity should occur at the position projected onto the 
star. Even though this is not the case here, the difference of the 
expansion velocity between the central slit and the Slit 8N is 
not significant and within the errors. 

However, a comparison with Slit 9S (PA = 50°) shows an 
asymmetric shape of the expansion structure indicating that the 
nebula is most likely not perfectly round but perhaps a tilted 
ellipsoid with the maximum expansion off-centered to the east 
of WRA 751. There the largest expansion velocities are found. 
Another possibility is that the nebula is spherical with a bump 
at this point. This is best visible in the asymmetric shapes of the 
expansion ellipses in Slits 8N (Fig. |) and Slit 9S (PA = 50°; 
Fig. While the largest negative velocity is found at the 0" 
position, the largest positive velocity is located at —2", i.e., the 
Doppler ellipse is asymmetric with respect to the projected po- 
sition of the star. All derived global expansion velocities are 
in good agreement with expansion velocities found in the lit- 
erature (Hutsemekers & van Drom 1991a, Garcia-Lario et al. 
1998). 
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In almost all spectra the redshifted line is more intense and 
broader than the blueshifted side. Since it is the redshifted com- 
ponent which is stronger, this brightness difference cannot be 
accounted by an absorption of the blueshifted component. It is 
either an intrinsic brightness difference or results from a thicker 
shell at the backside which leads to a higher emission measure. 
The FWHM of the redshifted wing of the expansion ellipse is 
resolved at 21 km s _1 (instrumental FWHM ~ 10 km s _1 ) and 
supports the interpretation as a thicker shell. The widths itself is 
indicative of a stratification of the radial velocities along a line 
of sight, presumably with the inner parts of the nebula moving 
slower. A wider redshifted component can be explained as due 
to a thicker shell there. 

Fig. H shows a composite pv-diagram of all slits with PA = 
140°. It shows the trend of the expansion velocities across the 
nebula: 

- The split of the Doppler ellipse decreases as we move 
away from the geometric center of the bubble. The de- 
crease, however, is not in agreement with a purely spher- 
ical expansion. While the redshifted part of the expansion 
ellipse decreases more like that of a spherical expansion, 
the blueshifted component stays nearly at a constant veloc- 
ity for a given position along the slit. 

- For a spherical expansion not only the width of the expan- 
sion ellipse in velocity space shrinks, but it also becomes 
narrower in spatial extent. In the WRA751 nebula, how- 
ever, we find that the convergence points at positive posi- 
tions migrate to smaller values as one proceeds to slits far- 
ther away from the star while they are at approximately the 
same location for negative values. 

Both points are in agreement with an asymmetry within the 
nebula and indicate deviations from a spherical expansion. 

Besides the — asymmetric — expansion of the shell, the 
spectra show additional kinematic components. At slit position 
8N the redshifted side converges with the blueshifted side in 
the north-west, at the position of about 10" as expected at this 
position for an expanding bubble. At the same point, however, a 
component which is more redshifted as the nebula appears, and 
extends to a position of 13". The redshift of this component 
increases to 32 km s _1 above the central expansion velocity as 
one moves away from the nebula. In Fig. || a dashed line in 
the ^-diagrams for Slits 8N, center, and 8S indicates the rest 
velocity of the kinematic center of the expansion at lOkms -1 , 
14kms _1 , and 12kms _1 , respectively. 

A similar kinematic behavior can be found in Slits cen- 
tral and 8S, but this time at the south-eastern rim of the neb- 
ula. At about the position —12" in both spectra an additional 
blueshifted component appears with maximum velocities of 12 
and 18kms _1 below the respective central expansion veloci- 
ties. Again the velocities of these components reach their max- 
imum at the location most distant from the star. 

The position of the two kinematic extensions is very sym- 
metric with respect to the central star and thus indicative of a 
bipolar kinematic structure. At Slit 8N material at the rim of 
the inner (nearly) spherical nebula moves away from us while 
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Fig. 5. Overlay of all position velocity diagrams at PA = 140°. 

in Slit 8S material at the corresponding positions, namely the 
edge of the nebula, approaches us. With respect to the center 
of expansion of the central spherically expanding nebula these 
kinematic components are bipolar. 

The spatial distribution of the velocity field corresponds 
well to the morphology discussed above: The main part of the 
expansion (between the two convergence points of the expan- 
sion ellipse) measures 23" in diameter and thus agrees with the 
diameter of the main body of the nebula as determined from the 
image to be 22" 8. Moreover, it is noteworthy that the sizes of 
the kinematic extensions as measured from the spectra and of 
the Northern and Southern Caps as determined from the images 
agree equally well. 

To ensure that the Caps are part of the nebula around 
WRA751 and not just a background object onto which 
the nebula is projected, we have determined the line ratio 
[Nil] 6583A/H a after subtracting contributions of the telluric 
H a line. Within the accuracy of our measurements, the en- 
tire nebula, including both Caps, shows a constant value of 
[N II] 6583A/H a = 1.2 ± 0.1, in contrast with the value for the 
background material of 0.3. The background line ratio is com- 
patible with the value for a galactic H II region (Shaver et al. 
1983). This difference ensures that we can disentangle source 
and background contributions. 

5. Discussion and conclusions 

A comparison of detailed kinematic data and the morphol- 
ogy of the LB VN around the galactic star WRA 75 1 indicates 
that the expansion is not as perfectly spherically symmetric as 
previously thought. The redshifted central part of the nebula 
expands nearly spherically while the blueshifted part appears 
stalled and moves at an almost constant radial velocity of about 
— 15 km s _1 (Fig. ^]). A natural reason may be that the front 
side of the nebula expands into a denser medium and that thus 
the expansion is decelerated. The maximum expansion veloc- 
ity was found off-centered from the central star adding a further 
piece of evidence that the nebula is neither exactly spherical nor 
expanding accordingly. The off-centered maximum velocity is 
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Fig. 3. Echellograms (left column) and corresponding position- velocity diagrams (right column) for the slits with PA = 140°. 
The positions are centered on the projected location of the star, the velocities are with respect to the LSR. Dashed lines in the 
pv-diagrams mark the center of expansion velocities. 
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Fig. 4. Echellograms (left column) and corresponding position- velocity diagrams (right column) for the slits with PA = 50°. The 
positions are centered on the projected location of the star, the velocities are with respect to the LSR. 



either an evidence for an elongated (elliptical) shape and tilt or 
for a bump on the back side of the nebula. 

Two kinematic extensions, one redshifted, one blueshifted 
appear attached to the ridge of the nebula. A comparison with 
the morphology of the nebula around WRA 751 allows to iden- 
tify the redshifted extension (Slit 8N) with the Northern Cap 
and the blueshifted one (Slits center and 8S) with the Southern 
Cap (see Cap and Slit positions in Figs, [j] and||). The North- 
ern and the Southern Cap therefore indicate morphologically 
and kinematically bipolar components in the nebula around 
WRA 751. While one may speculate about a jet like structure 
or — similar to planetary nebulae (for a discussion of similar- 
ities between LBVs and planetary nebulae, see, e.g., Frank 
1998) — wider funnels, or an altogether elliptical form of the 
nebula creating the bipolar appearance, due to the limited spa- 
tial resolutions of both the available kinematic data and the 
available images, higher resolution observations will be nec- 
essary to further constrain the structure of the nebula around 
WRA 751. 

Nonetheless, from the data presented one can draw several 
conclusions: 

- The LBVN around WRA 751 consists of an expanding 
shell of 23 " in diameter, with a thicker shell at the back 
side. The radial velocity distribution deviates from that of 
a simple spherical expansion pattern. We find bipolar con- 
tributions to the morphology as well as radial velocity dis- 
tribution in approximately north- south direction (Northern 
and Southern Caps). 



- The maximum expansion velocities are off-centered which 
most likely can be accounted for by a bump in the back side 
east of the central star. 

- The [Nil] 6583A/HQ, ratio is considerably larger than that 
of the background material. This is typical for CNO- 
processed material. 

These three results are fully consistent with WRA 751 being 
a true member of the LBV class of stars. In many of the well 
investigated ones bipolar components to the morphological as 
well as kinematic structures of the nebula are found, r] Car 
(Duschl et al. 1995, Morse et al. 1998) and HR Car (Weis 
et al. 1997, Nota et al. 1997) are the most prominent exam- 
ples. As stars in late phases of their evolution, LBVs quite 
naturally show CNO-processed material in their envelopes and 
ejecta (Garcia-Segura et al. 1996, Smith et al. 1978). A large 
[N II] 6583A/HQ, ratio in a nebula is often used as one of the in- 
dicators for a star being a LBV. One finds, for instance, values 
of [Nil] 6583A/H a = 3 ... 7 for rj Car (Davidson et al. 1982, 
Meaburnet al. 1987, 1996, and Weis et al. 1999), = 0.4 . . . 0.9 
for HR Car (Hutsemekers & van Drom 1991b, and Weis et al. 
1997), and = 0.7 for AG Car (Thackeray 1977 and Smith et al. 
1997). 

Combining our morphological, kinematic and spectro- 
scopic results, we find mounting evidence that WRA 751 is 
a LBV, indeed. Moreover, there are bipolar components in its 
nebula, notably the two Caps. While the bipolarity in this ob- 
ject is less pronounced and less obvious than in other LBVs, it 
still strengthens our suspicion that bipolarity — albeit at differ- 
ent levels — is a genuine property of LBV nebulae. This makes 
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it even more likely that modeling LB VNs as simple windblown 
spheres is an unjustified oversimplification. 
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